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I. INTRODUCTION 
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Formalism based on complex-scaling method is developed for solving the few particle scattering 
problem by employing only trivial boundary conditions. Several applications are presented proving 
efficiency of the method in describing elastic and three-body break-up reactions for Hamiltonians 
which may include both short and long-range interaction. 
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■ The quantum-mechanical problem of interacting particles is of the fundamental importance in theoretical physics, 
opening broad field of application related with the description of the microworld. Quick development of the computa- 
tional techniques, following the rapid evolution of the computational power, provoked sizeable advance in multi-particle 
bound state problem: rigorous and accurate description of the systems formed of up to several or even dozens of par- 
ticles have been obtained [J-Q- The progress in multi-particle scattering problem is moderated however. The major 
obstacle turns to be the rich variety of the reactions one should consider simultaneously and the resultant complexity 
of the wave function asymptotic structure. Till now only three-body system has been treated in a full extent, com- 
prising elastic and break-up channels 043; whereas rigorous description by the same methods of the four particle 
scattering remains limited to the elastic and rearran gem ent channels |8l-[l0| . Recently very courageous effort has been 
undertaken to apply Green Function Monte Carlo |lTI| and No-Core Shell model calculations for the nucleon 
' scattering on A>4 nuclei, nevertheless these promising techniques remains limited to the description of the binary 
scattering process. Therefore a method which enable the scattering problem to be solved without explicit use of the 

^ asymptotic form of the wave function is of great importance. 

I The complex scaling method has been proposed [13, [3 ^'^'^ successfully applied for the resonance scattering [l^ , 
CsJ , as has been demonstrated recently this method can be extended also for the scattering problem [l^ |T3l . In this study 
' we propose novel method to solve quantum few-particle scattering problem based on complex scaling method, which 
, allows to use trivial boundary conditions. We demonstrate success of this method in both calculating elastic and 
three particle break-up observables. 



II. FORMALISM: 2-BODY CASE 

1. Short range interaction 

The complex scaling method has been proposed a while ago to treat the scattering problem for the exponentially 
bound potentials [13]. Idea is quite simple and can be summarized as follows. First one recast the Schrodinger 
equation into inhomogeneous (driven) form by splitting systems wave function into the sum 4* = 4''"^ + containing 
the incident (free) \E'™(r) = exp(ik • r) and the scattered vl'*'^(r) waves as: 

{E-Ho- V{r))^"'{r) = y(r)*''"(r). (1) 

The scattered wave in the asymptote is represented by the outgoing wave ^'^'^oo exp{ikr) jr. If one scale all the particle 
coordinates by a constant complex factor, i.e. Yi = e^^Vi with Im{e'^^) > 0, the scattered wave vanish exponentially as 
^' ooexp(— fcrsin^) with increasing particle separation r. Moreover ff the interaction is of short range (exponentially 
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bound), the right hand side of eq. ([T]) also turns to zero for large r, enabling to solve former equation in a similar way 
as a bound state problem: using compact basis or by solving differential equation on a finite domain and by imposing 
to vanish on its borders. 
In practice one solves 2-body problem by expanding Schrodinger equation into partial waves: 



2/x 



Hoi{r)-Vi{r))rr{r)^Vi{r)i;r{r), 



(2) 



The radial part of the incoming wave is represented by the regular Bessel functions ipY^{r) = ji{kr)kr and with kinetic 
energy term given by 



!;2 r 



after the complex scaling, this equation becomes: 



2^ 



(3) 



(4) 



The complex scaled inhomogeneous term is easily obtained by using analytical expressions for the regular Bessel 
function tpi (r) = ji{kre^^)kre^^ . Extraction of the scattering phase-shift may be done either directly by determining 
asymptotic normalization coefhcient of the outgoing wave: 



(r) = Ai{r) exp{ikre - ln/2) 

with scattering amplitude given by: 

kfi = e'^' smdi= Mr ^oo) 
Other well known alternative is to use integral representation, which one gets after applying Green's theorem: 

f Jiikre'')Vi{re''')(^;\r)+i^Tir))re^''dr 



(5) 
(6) 

(7) 



2. Coulomb plus short range interaction 



If interaction contains long range term the problem turns to be quite different. The right hand side of m. ([T|) after 
the complex scaling diverges and the ^''^'^(re*^) term is not anymore exponentially bound. In ref. [l6l. [l7j exterior 
complex scaling was proposed as a solution to circumvent the problem due to diverging term on the right hand side 
of eq. (HJ. 

In this paper, when considering problem of the interaction containing short range {V'^{r)) plus Coulomb term 
{V^ (r) — ^j^) we propose to stick with the standard smooth scaling procedure, however to employ analytically 
continued Coulomb waves to circumvent the problem of the divergence. In this case driven partial wave Schrodinger 
equation is written as: 

- Hoi{r) - V%r) - {r))^^'' (r) = F^OC^^M, (8) 

where i/jj"''^ — Fi{r],kr) is a regular Coulomb function, which is solution of the former Hamiltonian containing 
Coulomb interaction only. Asymptotically the scattered wave behaves as ip'''^'^ oo exp{ikr — ri\n2kr), and therefore 
vanish exponentially after the complex scaling. 

^f^(r) = Aiir) exp{ikre'^ - 'qln2kre''^ - Itt/2). (9) 



Equation (|8]) may be readily solved with the vanishing boundary condition for 'ip^^{r), one simply must be able to 
continue analytically the regular Coulomb functions standing on the right hand side. 

The scattering amplitude and Coulomb corrected phase shift due to short range interaction 5i can be determined 
as previously from the asymptotic normalization coefficient 
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""7i = e^(*'+'"') sinSi = Ai{r oo), (10) 



where ai is so-called Coulomb phase shift. 

Alternatively, Green's theorem may be used to obtain integral relation similar to eq.©: 



fi^-^e''"' I Fi{r],kre'")V''{re"')^ {rye^'Ur (11) 

III. FORMALISM: 3-BODY CASE 

3. Short range interaction 

In sake of simplicity let consider a system of three identical spinless particles submitted to short range pairwise 
interactions. Only two vector variables are needed in the barycentric system, which may be one of the Jacobi pairs 
Xi = Tj—Tk and = (r^ — (r^ + r^)). The pair potential is assumed to support any number of two-particle bound 
states 0^(xi) with eigenvalues £,« and the angular momentum of this state Im- The corresponding continuum state 

has relative momenta q„i, satisfying energy conservation relation E = + = K^', the second equality defines 

three-particle break-up momenta K. 

Three particle problem we formulate by using Faddeev equations [31 in configuration space and by readily separating 
incoming wave of the particle scattered on a bound pair in the state cfy^^^'x.i). Three- identical particle problem maybe 
concluded in a single equation: 

[E-H^- y.(x,))V- ™(x., y,) - V,{^,) rj%,i^j,ys) = ^.(xO J2 '^r. (x,) exp(*q„, • y,), (12) 

where V'l^m is the scattered part of the Faddeev amplitude, corresponding to the incoming particle i; while by Vi we 
denote a pair interaction of the particles j and k. Decomposition of the systems wave function into three Faddeev 
amplitudes permits to separate two-cluster particle channels, whereas three-body break-up component remains shared 
by the three Faddeev amplitudes. In the yi — >■ oo asymptote the scattered part of the Faddeev amplitude i takes form 

V'i,m(Xj,yj) = Am{x^,y^,Xi/yi) ^^^^!^'^^ + V /m« (yi)0„ (x» ) '^^^ , '^"^'-^ (13) 



here p = y/xf+yf is hyperradius of the system; Am(xi,yi,Xi/yi) is a three-particle break-up amplitude; fmniyi) 
is two-body transition amplitude from channel m to channel n. In this expression sum runs over all open binary 
channels n. 

One may easily see that the scattered part of the Faddeev amplitude might vanish for large hyperradius if particle 
coordinates are properly complex scaled and thus: Xi = Xie*^,yj = y^e'^ and p = pe'^. However in order to perform 
resolution of the problem on a finite grid one should ensure that the inhomogeneous term, standing in the right hand 
side of eq. (jl2p . also vanishes outside the resolution domain. The inhomogeneous term is null, damped by the potential 
term, if Xi is large and thus falls outside of the interaction region. Alternatively for Xi <C j/i, the modulus of the 
transformed Jacobi coordinates approach Xj ~ ^yi] yj ~ yi/2 and 



,9 exp(-jg„yje*'') \/3 exp(9„^sm^ 
T-r^F oc exp -fc„— y,cos6l) ^ 



</>„(x,e^«)exp(zq„ . y,e^% oc (x,e'^) ""^^^^j^^f/'^ ' oc exp(-fc„^y, cos 0)^^il^:^li^^, (14) 



here we explored the fact that the bound state wave function decrease exponentially in the asymptote with momenta 
kn = \/ TTi \ tn\l'h. The last expression vanish for large yi only if condition 



Qn \l E + |e„| 

is satisfied [1^ . Therefore if the scattering energy is large enough one is obliged to restrict the complex scaling angle 
to very small values 
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Extraction of the scattering observables maybe realized as previously in two different ways. Straightforward way 
is to extract transition amplitudes from the — ;> oo asymptote of the solution tp^ (xi,yi), exploiting the fact that 
different scattering channels are mutually orthogonal: 

/„„(y,) = liMexpi-tq^y.e'') J (x,e"''')^;;„(x„ y^e^'^d^x,, (16) 

where C„ is normalization coefficient of two-body wave function f 0* (xie~*^)0„ (xie*^)e'^*^(i^Xi = C„. Break-up 
amplitude might be extracted from the tp^ (x.^, y.^) once all the two-body transition amplitudes are calculated relaying 
on eq. p6)) . 

Alternatively one can employ Green's theorem. In this case integral relations might be obtained both for break-up 
and two-body transition amplitudes. For the transition amplitude one has: 

/_(^,) = -C-'^ J J 0:(x,e-«)^^^P^^^p^(T/,(x,e^'^) + Ffc(x,e^«))W„,(x„ y^e^'^rf^x.dV. (17) 
These integrals are convergent on the finite domain, if the following condition is satisfied: 



--^^ - ^^"^ - , (18) 
qm + 2qn ^\e^\+E + 2y/\en\+E 

which is stronger than condition eq. (jlSp 26] . 

For the break-up amplitude several different relations can be obtained [l9| , it seems that the one employing 2-body 

outgoing states (/>*'~'"^(p, Xj), generated by the correspondingly scaled strong potential at relative momenta p, seems 
the most reliable numerically: 



In practice, calculations are performed by expanding former equations into partial waves. This pure technical issue 
is not subject of this paper and we refer interested reader to (l9| for the details on the partial wave Faddeev equations. 
One should note that partial wave expansion has no effect on the validity of the presented method. 



^. Short range interaction plus Coulomb 



The former discussion can be readily extended for the case when particles interact via short range plus Coulomb 
forces. In this case we prefer to use Faddeev-Merkuriev equations [1^, which elaborate Faddeev formalism by accom- 
modating Coulomb force. Indeed, Faddeev equations, suppose free asymptotic behavior of the particles, in case of 
long range interaction become ill behaved due to noncompactness of their kernel. These equations can still provide 
satisfactory solution for bound state problem, but are impractical for the scattering case [13] . Faddeev equations does 
not shed light on the asymtotic behavior of the separate amplitudes when long range interaction is present. Idea of 
Merkuriev [20j is to split the Coulomb potential V'-^ into two parts (short and long range), = V^''^ + V'"'-' , by 
means of some arbitrary cut-off function x- 

Vf-''{x,,y,)^Vf{x,)x^{x,,y,) V^-^ ix,,y,) ^ Vf" {x,)[l ~ x^ix^,y^)] (20) 
and reshuffle long range terms. One is then left with a system of equivalent equations: 

^E-Ho-W.- F/)V^r - Vf ^(V'f + C™) = (^^ - ^1 - Vf-^'m^l'^ = + Vj + Vt (21) 

by V^ we consider sum of the short range interaction Vi{xi) of the pair (jfc) plus short range part of the Coulomb 
force Vf = Vf''-^ + Vi. The incoming state is defined by: 

V'™„(x„y,) = 0„(x,)^^((?™,y,) (22) 

where (p^ {qm,yi) is a plane wave of the incident particle i moderated by its residual Coulomb interaction with a 
cluster of particles (jfc): Vf''''^'^''{yi) — -^j^^- The 0,^^(xi) is the eigenfunction of the m— th bound state of the 
particle pair (jk). 
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Former equations are solved as well as the scattering observables extracted in a similar way as for the Coulomb free 
case. For example transition amplitude, via Green's theorem, is expressed as: 

with Vi{xi) representing full interaction between particle-pair (jk). 

There is however formal difficulty associated with the extraction of the break-up amplitude for the case when all 
three particles are charged, since the asymptotic form of the break-up wave function is not known. One may still rely 
on the approximate relation employing Peterkop integral [21| as is claimed in [22i] . 

IV. RESULTS 

To test our approach we consider model of nucleons with mass ^ = 41.47 MeV ■ fm?^ where strong part of 
nucleon-nucleon (NN) interaction is described by S-wave MT I-III potential, defined as: 

Vs{r) = (Asexp(-1.55r) -f 1438.72 exp(-3.11r))/r (24) 

where Vs{r) is in MeV and r is in fm units. The attractive Yukawa term is defined with As=o = —513.968 MeV ■ fm 
and As=i — —626.885 MeV ■ fm for the two-nucleon interaction in spin singlet and triplet states respectively. 

MT I-III potential has been chosen for two reasons. On one hand it is widely employed potential for which accurate 
benchmark calculations exist. On the other hand this potential, being a combination of the attractive and repulsive 
Yukawa terms, reflects well the structure of the realistic nucleon-nucleon interaction: it is strongly repulsive at the 
origin, however has narrow attractive well situated at r w 1 fm. Note that many numerical techniques fail for the 
potentials, like MT I-III, which contain a repulsive core. 

In figure[T]we present our calculations for the NN ^Sq phaseshift at Ecm MeV. Two calculation sequences have 
been performed by enforcing the tpi ('') to vanish at the border of the numerical grid set at 50 fm and at 100 fm 
respectively, whereas complex scaling angle {9} has been chosen to be 10° and 30°. The phaseshift is extracted from 
the ipi (r) value at fixed distance r, according to eq. [5]|6] (Coulomb free case) and eq. [MTOl (short range plus Coulomb 
interaction) . Extracted phaseshift oscillates with r - this oscillatory behavior is due to the " premature" enforcement 
of ipi (r) to vanish at the border of the grid rmax- The phaseshifts extracted close to rmax are strongly affected by the 
cut-off and thus not reliable. The amplitude of the close-border oscillations is sizeably reduced by either increasing 
''max or 9, i.e. by reducing the sharpness of the numerical cut-off. Extracted phaseshift from the calculation with 
^max = 100 fm and 9 =30° is stable in a rather large window, which starts at r ~ 5 fm (right outside the interaction 
region) and extends up to r ^ 70 fm (beyond this value effect due to cut-off sets in) . In the stability region extracted 
phaseshift value agrees well with the exact result (dotted line). 

In figure[2]we compare NN ^Sq pahaseshift calculation for Ecm—^, 5 and 50 MeV by setting r„iax = 100 fm and 9 
=10°. One may see that by increasing energy effect of the cut-off reduces, sizeably improving stability of the extracted 
phaseshift. Inclusion of the repulsive Coulomb term does not have any effect on the quality of the method. 

One may improve considerably accuracy of the phaseshift calculation by employing integral relation eq. ([7]), see 
tables mill and figure |3l Phaseshift converges to constant value by either increasing cut-off radius rmax or complex 
rotation angle. Accuracy of five digits is easily reached. One should notice however that the use of very large values 
of 9 should be avoided, due to the fact that the function ?/'; (r) as well as complex scaled potential V{re^^) might 
become very steep and rapidly oscillating. At higher energy the function ipl'^ {r) vanishes faster and thus one achieves 
convergence by employing smaller values of rmax and/or 9. 

Our analysis has been extended to the three-body case. We consider nucleon-deuteron (N-d) i = scattering 
in spin-doublet [S = 1/2) and quartet {S = 3/2) states. Calculations has been performed both below and above 
three-particle break-up threshold. Below the break-up threshold results are stable and independent of the scaling 
angle, similar to the two-body case. Phaseshifts might be accurately extracted both using differential and integral 
expressions. 

Application of the differential relations for extracting scattering phaseshift and inelasticity above the break-up 
threshold is not so-obvious. It is difficult to find the stability domain. Therefore we employed integral expressions 
eqs. ([23l) and (fT9|). obtained using Greens theorem, which once again proved their worth. We summarize obtained 
results in Table IIIII and IIVI respectively for n-d and p-d scattering above the break-up threshold. Very accurate 
results are obtained for both phaseshift and inelasticity parameter once complex scaling angle is chosen in the interval 
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TABLE I: Calculation of the scattering phaseshift using integral expressions at Ecm =1 MeV 





MT I-III 


MT I-III+Coulomb 


rmax (fm) 


5° 


10° 


30° 


50° 


5° 


10° 


30° 


50° 


10 


44.420 


49.486 


55.790 


56.676 


33.999 


36.390 


41.528 


43.805 


25 


34.704 


44.211 


62.654 


63.743 


24.772 


34.910 


50.693 


50.698 


50 


56.812 


61.083 


63.482 


63.512 


39.895 


46.546 


50.487 


50.491 


100 


66.502 


63.822 


63.512 


63.512 


55.463 


50.811 


50.491 


50.491 


150 


62.497 


63.485 


63.512 


63.512 


49.317 


50.474 


50.491 


50.491 


exact 


63.512 


50.491 



TABLE II: Calculation of the scattering phaseshift using integral expressions at Ecm =50 MeV 



Tmax (fm) 


MT Mil 


MT Llll+Coulomb 




3° 


5° 


10° 


30° 


3° 


5° 


10° 


30° 


10 


19.400 


19.719 


19.923 


19.605 


19.795 


20.245 


20.610 


20.313 


25 


20.788 


20.135 


20.027 


20.032 


21.530 


20.864 


20.755 


20.760 


50 


20.014 


20.026 


20.027 


20.027 


20.734 


20.754 


20.755 


20.755 


100 


20.027 


20.027 


20.027 


20.027 


20.755 


20.755 


20.755 


20.755 


exact 


20.027 


20.755 




FIG. 1: (Color online) ^So NN phaseshift at Ecm~l MeV extracted using relations eg. I5l6l and [9lllOI respectively. Calculations 
performed with cut-off imposed at rmax=50 and 100 fm using complex rotation by the angle 8 —10° (dashed lines) and —30° 
(solid line). Pure strong interaction result is presented in the left figure, calculations including repulsive Coulomb interaction 
for pp-pair are presented in the right figure. 



[4°, 12.5°] for incident neutron with energy EiQb=14.1 MeV and [3°, 7.5°] at E;ab=42 MeV. Stability of the final resuh 
within at least three significant digits is assured, providing excellent agreement with the benchmark calculations 
of [IM [131 ■ Calculated integral gradually cease to converge on the finite domain for the calculations when higher 
complex scaling angles are chosen. This happens due to the fact that condition eq. p^ . which set limit 6'„iaa;=:14.2° 
and 9,nax=8.9° for the calculations at Eiaf,=14.1 MeV and Eiob=42 MeV respectively, is violated. 

In the table |V] we tabulate '^S'l n-d break-up amplitude as a function of the break-up angle -d, which defines pair 
and spectator wave numbers via k=K cos(t?) and q= 2K sin('i9)/-\/3 respectively. Nice agreement is obtained with the 
benchmark calculation of |23|] . Small discrepancy appears only for the i? values close to 90° , which defines configuration 
where one pair of particles after the break-up remains at rest. This is due to the slow convergence of the integral 
ea.(|19p for "i? — 90° in y-direction, special procedure must be undertaken in this particular case to evaluate the part 
of the slowly convergent integral outside the resolution domain defined by jmax ■ 
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FIG. 2: (Color online) ^So NN phaseshift calculation at Ecm~l, 5 and 50 MeV extracted using relations eg. I5I6I and [QlflOl 
respectively. Calculations performed with cut-off imposed at rmax= 100 fm using complex rotation by the angle =10°. Pure 
strong interaction result is presented in the left figure, calculations including repulsive Coulomb interaction for pp-pair are 
presented in the right figure. 
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FIG. 3: (Color online) Dependence of the calculated NN ^So phaseshift using integral expression as a function of the complex 
rotation angle. Grid was limited to rmaa: = 100 fm. The upper curve correspond Coulomb-free case, the bottom one includes 
Coulomb. 



V. CONCLUSION 



In this work we have presented a method based on the complex scaling, which enables to solve few-nucleon scattering 
problem without explicit treatment of the boundary conditions using square-integrable functions. Validity of the 
method is demonstrated for two and three particle scattering, including the three-particle break-up case with repulsive 
Coulomb interaction. Three-digit accuracy maybe easily obtained using this method. 

As is well known for two body case complex scaling angle is, in principle, only limited to 90°. On the contrary in 
order to solve three-body break-up problem the scaling angle should be restricted stronger from above, see eg. dTSl) . 
I.e. for the scattering at high energy the scaling angle should be limited to very small values. Nevertheless this 
limitation does not spoil the method at high energies, since rapid vanishing of the outgoing wave after the scaling is 
ensured by the large wave number values. 
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TABLE III: Neutron-deuteron scattering phaseshift and inelasticity parameter as a function of the complex rotation angle 6 
compared with benchmark results of [23. [24l|. Our calculations has been performed by setting ymai^lOO fm. 





3° 


4° 


5° 6° 


7.5° 


10° 


12.5° 


Ref. [23. 24] 








nd doublet at E;a6 = 


= 14.1 MeV 








Re{5) 


105.00 


105.43 


105.50 105.50 


105.50 


105.49 


105.48 


105.49 


V 


0.4559 


0.4638 


0.4653 0.4654 


0.4653 


0.4650 


0.4649 


0.4649 








nd doublet at Eiab'- 


=42 MeV 








Re(5) 


41.71 


41.63 


41.55 41.51 


41.45 


41.04 




41.35 


V 


0.5017 


0.5015 


0.5014 0.5014 


0.5015 


0.5048 




0.5022 








nd quartet at Eiab~ 


:14.1 MeV 








Re{S) 


68.47 


68.90 


68.97 68.97 


68.97 


68.97 


68.97 


68.95 


V 


0.9661 


0.9762 


0.9782 0.9784 


0.9783 


0.9782 


0.9780 


0.9782 








nd quartet at Eiab'- 


=42 MeV 








Re((5) 


37.83 


37.80 


37.77 37.77 


37.74 


38.06 




37.71 


V 


0.9038 


0.9034 


0.9032 0.9030 


0.9029 


0.8980 




0.9033 



TABLE IV: Proton-deuteron scattering phaseshift and inelasticity parameter as a function of the complex rotation angle 8 
compared with benchmark values of [2J|. Our calculations has been performed by setting ymai:=150 fm. 





3° 


4° 


5° 6° 


7.5° 


10° 


12.5° 


Ref. [24] 








pd doublet at Eiat- 


=14.1 MeV 








Re(<5) 


108.46 


108.43 


108.43 108.43 


108.43 


108.43 


108.42 


108.41 [3] 


V 


0.5003 


0.4993 


0.4990 0.4988 


0.4986 


0.4984 


0.4981 


0.4983[1] 








pd doublet at Eiab' 


=42 MeV 








Re{S) 


43.98 


43.92 


43.87 43.82 


43.78 


44.83 




43.68[2] 


V 


0.5066 


0.5060 


0.5056 0.5054 


0.5052 


0.5488 




0.5056 








pd quartet at Eiab~ 


=14.1 MeV 








Re{S) 


72.70 


72.65 


72.65 72.64 


72.64 


72.63 


72.62 


72.60 


V 


0.9842 


0.9827 


0.9826 0.9826 


0.9826 


0.9828 


0.9829 


0.9795[1] 








pd quartet at Eiab'- 


=42 MeV 








Re{S) 


40.13 


40.11 


40.08 40.07 


40.05 


40.35 




39.96[1] 


V 


0.9052 


0.9044 


0.9039 0.9036 


0.9034 


0.9026 




0.9046 
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TABLE V: Neutron-deuteron ^Si break-up amplitude calculated at Eiab=42 MeV as a function of the break-up angle i9. 





0° 


10° 


20° 


30° 


40° 


50° 


60° 


70° 


80° 


90° 


This work Re(^Si) 
This work Im(^Si) 


1.49[-2] 
1.69[0] 


8.84[-4] 
1.74[0] 


-3.40[-2] 
1.87[0] 


3.33[-2] 
1.92[0] 


7.70[-2] 
1.80[0] 


2.52[-l] 
1.68[0] 


4.47[-l] 
1.70[0] 


6.47[-l] 
1.96[0] 


6.30[-l] 
2.23[0] 


-1.62[-1] 
3.17[0] 


Ref. [23] Re(-'Si) 
Ref. [23] Im(^Si) 


1.48[-2] 
1.69[0] 


9.22[-4] 
1.74[0] 


-3.21[-2] 
1.87[0] 


3.09[-2] 
1.92[0] 


7.70[-2] 
1.80[0] 


2.52[-l] 
1.67[0] 


4.51[-1] 
1.70[0] 


6.53[-l] 
1.95[0] 


6.93[-l] 
2.52[0] 


-1.05[-1] 
3.06[0] 



9 



[5] J. L. Friar, B. F. Gibson, G. Berthold, W. Glockle, T. Cornelius, H. Witala, J. Haidenbauer, Y. Koike, G. L. Payne, J. A. 

Tjon, et al., Phys. Rev. C 42, 1838 (1990). 
[6] A. Kievsky, S. Rosati, and M. Viviani, Phys. Rev. Lett. 82, 3759 (1999). 
[7] A. Deltuva, A. C. Fonseca, and P. U. Sauer, Phys. Rev. C 72, 054004 (2005). 
[8] M. Viviani, S. Rosati, and A. Kievsky, Phys. Rev. Lett. 81, 1580 (1998). 
[9] R. Lazauskas and J. Carbonell, Phys. Rev. C 70, 044002 (2004). 
[10] A. Deltuva and A. C. Fonseca, Phys. Rev. Lett. 98, 162502 (2007). 

[11] K. M. NoUett, S. C. Pieper, R. B. Wiringa, J. Carlson, and G. M. Hale, Phys. Rev. Lett. 99, 022502 (2007). 
[12] P. Navratil, R. Roth, and S. Quaglioni, Phys. Rev. C 82, 034609 (2010). 
[13] J. Nuttal and H. L. Cohen, Phys. Rev. 188, 1542 (1969). 

[14] E. Balslev and J. M. Combes, Commun. Math. Phys. 22, 280 (1971)^ 

[15] N. Moiseyev, Physics Reports 302, 212 (1998), ISSN 0370-1573, URL'http : //www . sciencedirect . com/science/articl e/B6TVP-3VF47 
[16] M. B. C. W. McCurdy and T. N. Rescigno, J. Phys. B: At. Mol. Opt. Phys. 373, R137 (2004). 

[17] M. V. Volkov, N. Blan der, E. Yarevsky, and S_L. Yakovlev, EPL (Europhysics Letters) 85, 30001 (2009), URL 

|http : //stacks . iop . org /0295-5075/85/i=3/a=3 0001 [ 
[18] L. D. Faddeev, Sov. Phys. JETP 12, 1014 (1961). 
[19] W. Glockle, The quantum mechanical few-body problem (1983). 
[20] S. P. Merkuriev, Ann. Phys. (N.Y.) 130, 395 (1980). 
[21] R. K. Peterkop, Opt. Spektrosk. 153, 13 (1962). 

[22] P. L. Bartlett and A. T. Stelbovics, Phys. Rev. Lett. 93, 233201 (2004). 

[23] J. L. Friar, G. L. Payne, W. Glockle, D. Hiiber, and H. Witala, Phys. Rev. C 51, 2356 (1995). 

[24] A. Dehuva, A. C. Fonseca, A. Kievsky, S. Rosati, P. U. Sauer, and M. Viviani, Phys. Rev. C 71, 064003 (2005). 

[25] One must note that this condition is derived for a system of three particles with identical masses. For the case of three 

particles with arbitrary masses, one gets: tane < ^^J^^S^EIZI. where TTii is the mass of incoming particle, the binding 

energy En I correspond the weakest bound state of the particle pair {jk) open for scattering, M = rrn + mj + mk is a total 

mass of the system. 

[26] For the case of three particles with arbitrary masses, one gets: tan 6 < s/rriiM \e,n\/ {y^mjmk{\eni\ + E) + 

^{M - mj)(M - mk){\eu\ + E)) 
[27] In principle, Faddeev equations maybe solved also for long-range scattering problem after the complex scaling, however 

only vanishing of the total wave function and not of the separate Faddeev amplitudes (ones work with) is assured in this 

case. 



